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Abstract: 

Explicitly time-dependent, nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova rem- 
nants (SNRs) has been employed to investigate the properties of SNR RX J1713. 7-3946. Observations 
of the nonthermal radio and X-ray emission spectra as well as earlier H.E.S.S. measurements of the very 
high energy 7-ray emission were used to constrain the astronomical and the particle acceleration pa- 
rameters of the system. The model assumes that the object was a core collapse supernova (SN) with a 
massive progenitor, has an age of « 1600 yr and is at a distance of ~ 1 kpc. It is shown that an efficient 
production of nuclear CRs, leading to strong shock modification and a large downstream magnetic field 
strength £?d ~ 100 /iG, can reproduce the observed synchrotron emission from radio to X-ray frequen- 
cies together with the 7-ray spectral characteristics as observed by the H.E.S.S. telescopes. Small-scale 
filamentary structures observed in nonthermal X-rays provide empirical confirmation for this field ampli- 
fication scenario which leads to a strong depression of the inverse Compton and Bremsstrahlung fluxes. 
The results are compared with the latest H.E.S.S. observations. 



Introduction 

RX J1713-3946 is an extended shell-type SNR in 
the Galactic plane that was discovered in X-rays 
with ROSAT [15]. It turns out that the observ- 
able X-ray emission is entirely nonthermal. The 
radio synchrotron emission is weak, with a poorly 
known spectral form. RX J1713-3946 was also 
detected in very high energy (VHE) 7-rays with 
the CANGAROO [14, 9] and H.E.S.S. [1, 2] tele- 
scopes. Especially the latter observations show a 
clear shell structure at TeV energies which corre- 
lates well with the X-ray contours from ASCA. 
The different observations are discussed in the pa- 
per of Berezhko & Volk ([7], hereafter referred to 
as BV06), on which we expand here. 

The basic aim of BV06 is an investigation of the 
acceleration of both electrons and protons in this 
remnant, the calculation of the nonthermal radi- 
ation spectrum of the source using an explicitly 
time-dependent nonlinear kinetic theory, and a the- 
oretical determination of the (spatial) morphology. 
In such a theory the particle acceleration process 



is coupled with the hydrodynamics of the thermal 
gas in the aftermath of the SN explosion. 

In comparison with other SNRs (SN 1006, Cas A, 
and Tycho's SNR) that were successfully described 
within the framework of this theory and finally al- 
lowed a prediction of the 7-ray flux (see [5] for a 
review), the case of RX J1713-3946 is more dif- 
ficult despite its spatially resolved 7-ray emission. 
First of all, astronomical parameters such as source 
distance, present expansion velocity and age are 
not well known. Still, following estimates from 
ASCA observations [12] and from NANTEN CO 
line measurements [10, 13] a rather close distance 
of d = 1 kpc appears indicated, and historical Chi- 
nese records suggest a low age of t = 1614 yr [17]. 
These parameters are consistent with the assump- 
tion of a core collapse SN of type Il/Ib which ex- 
ploded into the adiabatic, very diluted bubble cre- 
ated by the wind of a massive (M < 20 M ) 
progenitor star, whose likely compact remnant in 
the center is an identified neutron star (e.g. [8]). 
Secondly, the lack of knowledge of the spectral 
shape of the radio emission makes it difficult to 
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Figure 1: Spatially integrated spectral energy distributions of RX J1713-3946 [7]. The ATCA radio data, 
ASCA X-ray data, EGRET spectrum of 3EG J1714-3857, CANGAROO data and H.E.S.S. data from ([2]), 
are shown, together with the EGRET upper limit for the H.E.S.S. source position. The solid curve at 
energies above 10 7 eV corresponds to 7r°-decay 7-ray emission, whereas the dashed and dash-dotted curves 
indicate the inverse Compton (IC) and Nonthermal Bremsstrahlung (NB) emissions, respectively. 



derive - from synchrotron spectral observations 
- two determining physical quantities: the effec- 
tive strength of the magnetic field and the proton 
injection rate into the diffusive acceleration pro- 
cess. Strictly speaking, it is therefore not possi- 
ble to theoretically predict the TeV 7-ray emission. 
BV06 argue nevertheless that the observed overall 
synchrotron spectral shape, from radio frequencies 
to the X-ray cutoff, and the small-scale filamen- 
tary structures in the nonthermal X-ray emission 
of RX J 17 13-3946 are consistent with efficient CR 
acceleration associated with a considerably ampli- 
fied magnetic field. These properties allow in ad- 
dition a consistent fit of the observed TeV 7-ray 
spectrum. It is strongly dominated by 7r°-decay 
emission [2]. 



Model 

For a present angular radius of the SNR of 60', 
at the adopted distance of d = 1 kpc the lin- 
ear remnant radius corresponds to 10 pc. To pro- 
duce the high observed 7-ray flux the SN shock 
must already propagate into the increasing den- 
sity TVh of the shell of swept-up ambient inter- 
stellar medium (ISM) whereas the high nonther- 
mal X-ray emission in the energy range of sev- 
eral keV requires a presently still high shock speed 
V s k 1840 km/s, and therefore a very low bub- 
ble density Nf, « 10~ 2 cm~ 3 . This requires an 
ISM density N IS m = 300 cm~ 3 , i.e. a bubble 
in a molecular cloud. A consistent hydrodynamic 
explosion energy is E sn = 1.8 x 10 51 erg, with 
an ejected mass of M e j = 3.5 M® which give a 
swept-up mass M sw w 20A/ Q M e j, already at 
the present time. 
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Figure 2: The same theoretical spectra above 10 7 eV as in Fig.l, on an expanded scale, together with the 
EGRET upper limit for the H.E.S.S. source position and with the latest H.E.S.S. spectrum from [3], 



The effective magnetic field strength B inside the 
SNR, taken as approximately uniform [6], deter- 
mines the strength and form of the electron dis- 
tribution for given observed synchrotron emission 
strength, and thus also the inverse Compton emis- 
sion at very high 7-ray energies for a given radi- 
ation field like the CMB. Since the spatially inte- 
grated radio synchrotron spectrum is not known in 
detail, B is determined using the existence of nar- 
row filaments in hard X-rays, whose thickness we 
interpret as the result of synchrotron losses of the 
highest-energy electrons behind the SNR shock. 
Recent XMM observations by Hiraga et al. [11] 
show such a narrow filament from which B V06 de- 
rived a lower limit to the interior field strength of 
B = 65/iG. It is clearly amplified relative to the 
field strength in the ambient shell and indicates ef- 
fective acceleration of nuclear particles [4]. The 
actual value used in the model for the synchrotron 
emission is w 130/iG. 



Results and discussion 

The amplitude of the observed VHE 7-ray emis- 
sion can be fitted with a proton injection rate 
T) = 3 x 10~ 4 , quite plausible from theoreti- 
cal considerations. This is consistent with the 
high B-field strength which implies a predomi- 
nantly hadronic 7-ray emission, since the leptonic 
channel is strongly suppressed for the given syn- 
chrotron emission. The overall shock compres- 
sion ratio is a = 6.3, i.e. the shock is signifi- 
cantly modified by the accelerated CR protons. It 
turns out also that our spherically symmetric ac- 
celeration model overproduces nonthermal particle 
energy E c . The necessary deviation of the SNR 
from spherical symmetry requires a renormaliza- 
tion f re E c with f re w 0.2, as for a similar object 
like Cas A. 

Our resulting overall nonthermal spectrum [7] is 
given in Fig.l. In the VHE range it shows the data 
from CANGAROO [9], together will the H.E.S.S. 
data [2] available at the time. The same theo- 
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retical 7-ray spectral energy distribution is com- 
pared in Fig. 2 with the latest H.E.S.S. measure- 
ments [3] which have significantly increased statis- 
tics and energy coverage. Despite the fact that 
the theoretical model does not contain any parti- 
cle escape at the highest energies, the spectrum 
falls off more quickly with energy than the mea- 
surement. We trace this back to the fact that the 
model of BV06 has - in a conservative sense - 
used the present value of B as a constant value for 
B during SNR evolution, whereas B 2 /(Sn) pre- 
sumably scales oc NhV 2 [6, 16] or even oc NhV^ 
[4]. These latter scalings would lead to a somewhat 
higher proton cutoff energy at the present epoch. 

Also the theoretical X-ray and TeV 7-ray morphol- 
ogy can be studied. As expected from the filament 
observations, the projected radial X-ray profile is 
very narrow. The corresponding profile at 1 TeV is 
again very narrow, since in particular the gas den- 
sity has a large radial gradient. Smoothed to the 
resolution of H.E.S.S., the profile is consistent with 
the observed H.E.S.S. profile. 

A major result remains the hadronic dominance in 
the 7-ray emission spectrum. We expect that the 
GLAST instrument will confirm our corresponding 
theoretical prediction that the spatially integrated 
7-ray spectral energy density at 1 GeV is only a 
factor « 2.5 lower than at 1 TeV. 
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